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A DISCRETE HARDY’S UNCERTAINTY PRINCIPLE AND DISCRETE 

EVOLUTIONS 

AINGERU FERNANDEZ-BERTOLIN 


Abstract. In this paper we give a discrete version of Hardy’s uncertainty principle, by using 
complex variable arguments, as in the classical proof of Hardy’s principle. Moreover, we give 
an interpretation of this principle in terms of decaying solutions to the discrete Schrodinger 
and heat equations. 


1. Introduction 


In functional analysis, uncertainty principles are, in general, results that state that a function 
and its Fourier transform cannot decay too fast simultaneously. The most studied uncertainty 
principle, which comes back to Heisenberg, says that 

^(f \xf{x)\'^dx\ ([ \yf{x)\'^dx\ > [ \f{x)\^dx, 

« Vjr<j / KJS.d- J Js.d 

and, moreover, the equality is attained if and only if f{x) = for a > 0. In the 

discrete setting there are some versions of this inequality. For instance, in a previous paper [5], 
we studied an inequality discretizing the position and momentum operators (see also [T1I21I1U] for 
more references to this uncertainty principle). We saw how we can recover the classical minimizer 
(i.e. the Gaussian) from the minimizer of the discrete uncertainty principle, which is given in 
terms of modified Bessel functions 

= - f cos{ke)d9, kGZ. 

Jo 

In this paper, we are interested in giving a discrete version of Hardy’s uncertainty principle 
HIS] in one dimension 

|/(x)| < |/(C| < with a/3 < 1 ^ / = 0, 

and, in the case aP = 1 then f{x) = Thus, this uncertainty principle is telling us that 

a function and its Fourier transform cannot have both Gaussian decay for some coefficients a, /3 
The original proof of this principle is strongly based on complex variable arguments (Phragmen- 
Lindelof’s principle and Liouville’s theorem). Moreover, if we consider a solution to the one¬ 
dimensional free Schrodinger equation 

f dtu = idxxU, 

\ u{x,0) = ■uo(x). 
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we can write the solution as 




^ / X 


y/ii 

Basically, this says that the solution is the Fourier transform of the initial datum, so we can 
write Hardy’s uncertainty principle in terms of solutions to the Schrodinger equation. 


( 1 ) 


u{x, 0) = 0(e 




), u{x, 1) = 0) = 0(e 




), a/3 < 4 => It = 0, 


having a precise expression for uq in the case a/3 = 4. On the other hand, it can also be stated 
in terms of solutions to the heat equation, as follows: 


( 2 ) 




/2S^dxa, j g for some (5 < 2 => / = 0. 


Under these terms, there is a series of papers [5HH] and [3] where the authors prove Hardy’s 
uncertainty principle for perturbed Schrodinger and heat equations just by using real calculus 
arguments. In the case of the Schrodinger equation, they are able to get a result up to the 
endpoint case. However, for the heat equation their result does not cover the whole case. 

Here first we are going to see, using complex variable arguments, a version of Hardy’s un¬ 
certainty principle in a discrete setting, saying that if a complex function is controlled in some 
region of the complex plane by a function that will be closely related to the Gaussian, and its 
Fourier coefficients are controlled precisely by the minimizer we got in [3] (i.e. by the modihed 
Bessel function /^(x)), then in some cases we have that the function is zero, or we can give a 
precise expression for the function. This is not exactly stated as Hardy’s classical principle, since 
now, by using Fourier series we relate the sequence to a periodic function in an interval of the 
form [—7r/h, 7r//i], while in the continuous case the relation one gets using the Fourier transform 
is between functions in the real line. Due to this periodicity, just knowing the behavior of the 
function in that interval does not give us any information. 

Once we have this result, we want to give some uncertainty principles in terms of the free 
Schrodinger and the heat equations, concluding that solutions to these equations cannot decay 
too fast at two different times. We are going to use two different approaches to prove these 
results. First, we take advantage of what is known in the continuous case, so we need to assume 
that the initial datum of the discrete equation is similar to a function on the real line. In this 
case, it is not reasonable to expect that the solution to the discrete equation is identically zero, 
but it tends to the zero function in some sense as h, the mesh step, tends to zero. However, when 
this mesh step is fixed, we can use the discrete version of the uncertainty principle we mention 
above in order to see that under some cases this solution is identically zero. As it is pointed 
out in o, the condition on the decay coefficients in the continuous case is a/3 < 4, while, the 
discrete approach we are going to study here leads to a more restrictive condition for a and /3, 
a + /3 < 2. Although this is a clear difference between the continuous and the discrete case, we 
are going to see here that this is the sharp condition for our results. 

When preparing this manuscript, we learn about a recent and independent result in this 
direction m- There, the authors use complex variable arguments to prove a sharp analog of 
Hardy’s uncertainty principle considering solutions to the discrete Schrodinger equation. On the 
other hand, they also use the real variable approach in [SHE] in order to add real-valued time- 
dependent potentials. There is overlap between their results and those presented here in Section 
4, although the proofs of the results are different. 
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The paper is organized as follows: In Section 2 we give a general result in the discrete case that 
is related to Hardy’s uncertainty principle. In Section 3 we give a version of Hardy’s uncertainty 
principle for the discrete Schrodinger and discrete heat equations, by relating these equations 
to the continuous equations, noticing that we should not expect the solutions to be identically 
zero using this approach. In Section 4 we use the results of Section 2 to give discrete versions of 
Hardy’s uncertainty principle in the spirit of the results in Section 3 without using information 
about the continuous case and, moreover, we give some examples of discrete data that prove the 
sharpness of the decay conditions. Finally, in an appendix we give some examples of nonzero 
functions that satisfy the hypotheses of Theorem 2.1 when the statement of the theorem does 
not give any extra information about the function. 


2. A DISCRETE UNCERTAINTY PRINCIPLE USING COMPLEX VARIABLE ARGUMENTS 

For h > 0, that represents the mesh size, we define the Fourier coefficients of a ^-periodic 
function fh in the following way: 

A(fc) = ^ d^, fh{x) = ^ E 

‘SjI'K J — T^jh y I'K 


There is a result concerning the extension of this function jh to the complex plane. Roughly 
speaking, this result states that a function that decays faster than Va > 0, can be extended 

to the complex plane as an entire function. When the Fourier coefficients are of the form 1^ {u/h?) 
for u G C we are in this case, and we have that the corresponding periodic function is 


fcez 


^/^(y//j2)gi/tfc(a;+iy) _ _J]_^ucos{x+iy)/h 




So we are going to see that if the Fourier coefficients decay faster than the modified Bessel 
function, and on the other hand the periodic function fh is controlled by some function that 
is related to the one we have computed above, then in some cases the function will be zero, 
according to |u| and the argument of the Bessel function. 

Theorem 2.1. Assume that fh is a complex-valued function and that there are u = G C, b G 
[0, 27r), S G (0,7r/2), and s > 0 such that, for all y < 0, 

^ Ce^^ cos(^ — — (5) cosh^yh) — sin(0 — —5) sinh(yli) 

< Ce^ cos(S+.|+(5) cosh{yh)—f^ sin(e+.|+<5) sinh(i/?i) 


< Ce^ cos(0+.|+<5) cosh(y/t) + .^ sin(6'+.|+<5) sinh(i/?i) 
^ Ce^^ cos(^—(5) cosh(yIi) + ^^ sin(0——d) sinh(yh) 

Moreover, assume that the Fourier coefficients of fh satisfy 

\Mk)\<CIk(^-^y VfcGZ. 


fh 


fh 


b-e+^+s 


b-e-^-s 


and, for all y > 0, 
fh 


b 0 ^ S 


fh 


b+9-^-S 


ly 


ly 


ly 


ly 
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Then: 


(1) rs < 1 There are nonzero functions that satisfy the hypotheses. (See the appendix 
below) 

(2) rs = 1 ^ fh{z) = 

(3) rs > 1 ^ fhlz) = 0. 

Remark 2.1. The hypothesis on the periodic function is telling us that we know how the function 
behaves near the critical lines where the function e^ cos{zh-b) greatest decay. 

Remark 2.2. In the applications of this theorem, we are going to consider 6 close to 7r/2. 


Proof. To begin with, we show that the case rs > 1 can be proved as a consequence of the case 
rs = 1, as in Hardy’s original theorem. Indeed, if r > i, the hypotheses of the theorem are 

i6 

satisfied for m = ^ and s, and in this case rs = 1. To see this, let us consider y < 0, so that [C 
is a constant that may change from line to line) 


fh 


/b-0 + ^+S \ 

[ — h^^V 


<- cos(^— ^—S) cosh.{yh) — ^^ sin(0 — — 5) sinh(y/i) 

^ cos(7r/2+(5)/2/i^ ^ cos{7t/ 2-\-S)/2sh'^ 

<- cos(^— ^—8) cosh.(yh) — ^^ sin(0 — — 5) sinh(y/i) 


since y is negative and cos(7r/2 + (5) is also negative. For the other three lines that we have to 
control, we can use the same argument. Hence, by (2) in the theorem, we have that fh{z) = 

Ce'^ (( 2 -^')^) ^ hypotheses of the theorem are not satisfied, unless C = 0 ^ ft = 


0 . 


Thus, we only need to prove that when rs = 1 we can determine completely the function by 
the properties of the hypotheses. The proof of this fact relies on the maximum principle and 
Liouville’s theorem, so it is quite similar to the original proof of Hardy’s uncertainty principle. 
As we have pointed out above, the condition for the Fourier coefficients implies that fh is a 
27r//i-periodic function and the extension 

fh{z) =-^^fh(k)e^^'^\ z = x + iy, 
feGZ 

is an entire function, so we have, for all z in the plane 
\Mz)\ <ChY,h (^) 


We split up a strip of length 27r/h in six regions following the figures below (being the first 
one for the upper half plane and the second one for the lower half plane, and the boundary of 
each region is determined by the dashed red lines and the line = 0}): 
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0 


2 


3 


6+0—7r/2—(5 
h 


b±e 

h 


6+0+7r/2+^ 

h 


6+0+37r/2-(5 

h 




4 


5 




6 


6—0—7r/2 —5 
h 


b-9 


h 


6—0+7r/2+5 6—0+37r/2 —5 

h h 


The procedure is very simple. In each region, we first multiply fh by a nice function (which 
will depend on a parameter e) that decays when y tends to infinity in the unbounded part of 
the region. Then we see that the product is bounded at the boundary of the region, so we can 
apply the maximum principle and let e tend to zero. We do not want the sign of the real part of 
this function to change in each region, so we look for the behavior of sin a; cos a; on [7r/2,7r]. We 
illustrate this in region 1, since the procedure is exactly the same. In this region we consider the 
function 

g,{x + ly) = h{x + iy)e-i^ cos(./.-b)-..cos^(T.)^ 

where Tz = aiz + bi is the linear transformation which maps the interval into 

[7r/2,7r]. 


Hence, we have that 


Mx + iy) = + 


sin{xh — b) sinh{yh)-\-2€4>{x,y) 


where 4>{x, y) = cos(aia; + bi) sin(oia; + bi) cosh(aiy) sinh(ai?/). 


We can see that at the boundary of the region we are considering, this function is bounded. 
Moreover, a simple computation shows that + = > 1, and this fact tells us that when y 

tends to infinity, the leading term is given by the e part, which is negative. Hence the function 
is bounded when y is large and we can apply the maximum principle to get that | 5 e(z)| < C. 
Letting e tend to zero we conclude that, for all z in region 1, 


When we repeat this procedure in the other regions, we have that \fh{z)e~'^ < C 

in the two strips described in the figure above. One is a strip of length 27r/h of the upper half 
plane, and the other one is another strip of length 27:/h of the lower half plane. By periodicity, 
this estimate holds for all z in the plane, and by Liouville’s theorem this implies that f{z) = 

fj^-^cos{zh-b)_ g 
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In Section 4 we will use this theorem in order to see that a solution to a discrete equation 
cannot decay faster than some modified Bessel functions at two different times. What we have 
to see is the evolution of this equation in the periodic setting, by considering that the solution 
is given as the Fourier coefficients of some periodic function. The behavior at one time and the 
expression of the solution in the periodic setting will give us bounds for the function in some 
vertical lines of the complex plane, so we will use these bounds and the decay of the function at 
the other time to conclude that the function is zero. 


3. A DISCRETE Hardy’s theorem based on the classical Hardy’s uncertainty 

PRINCIPLE 

We are going to consider a solution to the Schrodinger equation 

( 3 ) + 

where, in order to use the continuous uncertainty principle in the discrete setting, we are going 
to assume that the initial datum /^(O) is not far from the evaluation of a continuous function 
uo{x), that we will use as initial datum for the continuous Schrodinger equation. More precisely, 
we require that there is a function uq such that 

Y,\fkiO)-uo{kh)\^ <h^, 
fcez 

for some fj. > 0. In order to be able to evaluate the function, we are going to require uq € 
for some s > 1/2. Moreover, this is telling us that at any tme t, the solution to the discrete 
Schrodinger equation with initial datum (0) is not far from the solution to the discrete equation 
with initial datum z^{0) = uo{kh), since the norm is preserved in time. 

Now what we want to see is how we can relate the solution to a discrete equation to a solution 
to its continuous version. Since we are going to give an £°° version of the uncertainty principle, 
we want to estimate the £°° norm of the difference between the evaluation of the solution to the 
continuous equation at time t in the mesh {kh : fc > 0} with initial data uo(x) and the solution 
to the discrete equation with initial data uo{kh). Since uq S with s > 1/2, we have that 

[ (l + |e|)(^^-'^/^|iio(e)Me<+oo, 

Jr 

and using similar arguments as in where this type of equations is studied when the time is 
also discrete, we can get that 

(4) \u{kh, 1) - z,^l)| < /■ (1 + dt 

Jr 

This is the easiest way to talk about convergence of the discrete equation to the continuous 
one. However, in m the authors give estimates for a vast variety of Banach spaces. 

Under these circumstances, we have the following result: 

Theorem 3.1. Let f^{t) he a solution to the discrete Schrodinger equation ([ 2 ]), and assume that 
there is uq G with s > 1/2 such that for some fi > 0, 

<h^- 

feez 


( 5 ) 
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If for some a, j3 satisfying aj3 < 4 we have, for all k and h > 0, 


( 6 ) 


l/fc^0)| < C 


Ik{a/h?) 


\fka)\<c 


hiP/hy 

lyp/hy-' 


then uo = 0 and \fjf{t)\ tends to zero uniformly in k when h tends to zero. 


Proof. Since we want to use the continuous result, what we need to prove is that uq and it (a:, 1) = 
are controlled by some Gaussians. We have to distinguish several cases according to x. 

First, if X = 0 we have that |/o (0) — ito(0)| < so 

\uom<\fUo)\ + h>^^^ <c+h>^/y 


and letting h tend to zero, we conclude that |ito(0)| < C. 

Now, if 0 < X < 1 we have that uo(x) = uq (jf ) for all j G N, and then Ifj^yO) — ito(a;)| < 
< - 7 ^- Hence 


|Mo(a:)| < \ ffy0)\ + 


ja /2 — lyaf^jxy 


1 

jM/2 ’ 


and it was proved in that this quotient of Bessel functions tend uniformly in x to the Gaussian 
g-rc / 2 q!^ gQ |;Q infinity we have |ito(a;)| < Ce~^ Changing j by —j we can 

argue in the same fashion to get the same result in — 1 < x < 0. 

Finally, if 1 < x we use that uo{x) = uq ^|~x] for all j G N, and then, as in the 
previous case 


|Mo(a;)| < 1/ 


\xYj I 


I\xY]{af\xY/xy 1 
|"x]^'^j^/^ — /o(aj^ |"x]®/x^) 


and we can use the same reasoning as in to see that this quotient also converges to the Gaussian 
uniformly in x. Actually, this case is easier than the previous one, since in that case one has to 
show a uniform estimate over compact sets |x| < M and then see that if M is large, for M < |x| 
both functions are very small. In this case, we do not need to do this and it can be proved the 
uniform convergence in one step. Now letting j tend to infinity we get that |ito(2;)| < Ce~^ 
in this region. Changing j by —j we prove the same when x < — 1 and gathering all the results 
we conclude that |uo(x)| < for all x G M. 


Now what we need to see is the evolution of uo(x) under the continuous Schrodinger equation. 
Since uq G i?®, with s > 1/2, we have that the solution to the continuous Schrodinger equation 
at the mesh is similar to the solution to the discrete Schrodinger equation when we discretize 
uo(x) by taking its values at the mesh, according to (H]). On the other hand, since the initial 
data /m(0) satisfies (O, that is, it is close to the evaluation at the mesh of uq, we have that at 
time t f^f{t) is also close to z!y{f), since we have 


\ft{f)-zt{t)\< 





1/2 



1/2 


Uo{kh)Y 


< 
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Hence, we can do the same as we have done with uo{x) to see that it is bounded by another 
Gaussian. For instance, if 0 < a; < 1, 


\u{x, 1)1 = |u ( j1 ) I < \Zj^\l)\ + 


j{2s-l)/8 


< i/;^'(i)i + 


j{2s-l)/8 ^ jf./2 - /^(y 3 j 2 / 3 . 2 ) 


1 

ji2s-l)/8 


1 

jp,/2 ’ 


and letting j tend to infinity we get |M(a;, 1)| < e ^ Doing the same for the other regions we 
finally get that \u{x, 1)| < e~^ for all x G R, so we can use Hardy’s uncertainty principle in 
its version for Schrodinger evolutions o to conclude that it = 0. In particular, this tells us that 
Mo = 0 and then what we have is that 1/^ (t)| < for all A: G Z, so the solution to the discrete 
equation tends to zero uniformly in k when h tends to zero. □ 


As we see in the statement of the theorem, we do not prove that the discrete solution is 
= 0, but it is close to zero as h is getting smaller. We can find examples where we see that 
we cannot expect to have a complete analog version to Hardy’s uncertainty principle using this 
approach. For instance, if we take as initial datum 


/fe (0) = 


Ik{il2h^) 

/o(5/2/i2) 


l/"(0)| < 


Ik{l/2h^) 

/o(5/2/i2) 


< 


4(l/2/i2) 

/o(l/2/i2) ’ 


since \Ik{z)\ < Iki\z\) and the modified Bessel function Io{x) is increasing when 0 < x < oo. 
Thus a = ^, and, if we solve the discrete Schrodinger equation we get at time t = 1, 


/fe(l)=e 




E 

m— — oo 


/o(5/2/i2) 


^k—m(2i/h ) — e 


2', _ „-2i/h' 


Imi5il2h^) 

/o(5/2/i2) 


l/fc (1)1 < 


/^(5/2/i2) 
/o(5/2/i2) ’ 


so in this case (i = \ and a/3 = | < 4. Then by the theorem, if there exists mq such that ([5]) 
is satisfied, uq has to be identically zero. In this case it is quite easy to see that this happens. 
Indeed, 


Ei/fc(o)i' 

fcGZ 


1 


and seeing the asymptotic behavior of Io{t) when t is large, we see that © is satisfied for any 
fj, > 0. Hence we have an example of a nonzero sequence that satisfies the hypotheses of the 
theorem. However, in the next section we will see that if we relax the condition on the coefficients 
a and /3, at least when h is fixed we can give a result that tells us when fj^ is going to be identically 
zero. 


On the other hand, if we consider as initial datum the sequence 


9k{0) 




then in this case we are not going to have a function mq satisfying since the subsequences 
of the odd members and even members are discrete version of the functions —and e~^ 
respectively, as we can see in Figure [TJ In this case, the solution to the Schrodinger equation is 


= £ (-1)’ 


(2ilhy - 


/o(l//l2) ’ 
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Figure 1. Graphic representation of g^{0) when h = -^ and k G [—50, 50] The 
odd coefficients are plotted in red dots and the even coefficients in black dots. We 
see the convergence of those subsequences to their respective Guassians ±e“^ 
represented in blue and green lines. 


and, at least numerically, the decay coefficient for g^{l) seems to be (see Figure^) P = 5. This 
is not very surprising, since the decay one obtains solving the continuous Schrodinger equation 
with initial datum is precisely , so the coefficients in this case are a = l,/3 = 5 as 

well. On the other hand, this is telling us that maybe one can remove the hypothesis ([S]) on the 
existence of a proper function uq and having then that just the decay conditions m imply that 
the solution to the discrete Schrodinger equation is getting smaller uniformly when h tends to 
zero. 


0.00003 f 
0.000025.1 

0.00002 - 

0.000015 - 

0.00001 - 

5.x 10“'’ - 


f f t 


Figure 2. Representation of 13^(1) | (black dots) when h = ^ and k G [200,250] 
We see that it seems to be controlled by for /3 = 5 (green dots), but 

this is not the case if we take /3 = 4.9 (red dots). 


Now, considering the discrete heat equation, we have similar results. Although in this case 
since we have a parabolic equation we can get better estimates for the convergence of the dis¬ 
crete system to the continuous equation, we are going to use the same estimates dl]) as for the 
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Schrodinger equation. Recall that we are comparing the discrete datum with the samples of 
a continuous function, and this is well-defined when the function belongs to H^,s > 1/2, so 
although it could be possible to give some estimates depending on the norm of the function in a 
weaker space for instance), the function has to have some degree of smoothness that makes 
useless the advantage of using these better estimates. On the other hand, in Theorem 3.1 we 
have used an estimate for the difference between the discrete initial datum and the continuous 
function, while now we can replace it by an i°° estimate. Indeed, if we assume that a sequence 
vl^ satisfies \vl^\ < for some ^ > 0 and for all m £ Z, then, solving the discrete heat equation 
we get 

CO CO 

k— — oo k— — oo 

so we have that the evolution of also satisfies the same estimate. In the case of the Schrodinger 
equation we do not have this property, and this is the reason why we have to consider an i'^ 
estimate, since we have the conservation of this norm. Then the result we have is: 

Theorem 3.2. Let v^{t) be a solution to the diserete heat equation, and assume that there is 
uq £ with s > 1/2 such that for some fi > 0, 

(7) sup\v^{0) — uo{kh)\ < . 

fcez 


( 8 ) 


If for some a satisfying a <2 we have, for all k £ Z and h > 0, 

Ikia/hf) 


kfc(l)l < C 




then uo = 0 and \v^{t)\ tends to zero uniformly in k when h tends to zero. 


Proof. Since we have that uq £ we only need to see if e® 1) belongs to L^(]R) for some 

5 <2. Since we have the uniform estimate © and the convergence of the solution to the discrete 
heat equation with initial datum z^ = uo{kh) to the solution to the continuous heat equation 
with initial datum uo{x) (|3]), we can repeat the procedure we use in Theorem 3.1 to conclude 
that \u{x, 1)1 < Now, since a < 2, que can take J such that a < 5 <2 and then 

f 1)1^ da: < / <-foo, 

Jr Jr 

so, by Hardy’s uncertainty principle for heat evolutions © we conclude that uq = 0. □ 


4. A WEAKER DISCRETE VERSION OF HARDY’S UNCERTAINTY PRINCIPLE 

As in the classical result of Hardy, we want to have decay conditions on solutions to the 
discrete Schrodinger equation that make the solution to be identically zero. For this purpose, 
we are going to use results coming from Section 2, where we use complex variable arguments to 
conclude that if a periodic function and its Fourier coefficients satisfy certain decay conditions, 
then the function is identically zero. 

The first question is to choose the parameters b and 6 in Theorem 2.1 properly. For that, 
we rewrite the discrete Schrodinger equation ([3]) in a periodic setting, thinking of (t) as the 
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Fourier coefficients of a function gh{t)- It is easy to see then that solving ([H]) is the same as 
solving 

dtgh{x,t) = => gh{x,t) = 0), 

and then we observe that solving the discrete Schrodinger equation is basically multiply by the 
exponential so we are going to take b, 9 so that the conclusion of Theorem 2.1 is that if 

rs = 1, then the function is a constant times e’^ cos(zft)^ that is, we take b — 0 and 6 = ^. Then 
we can rewrite Theorem 2.1 in order to have: 


Corollary 4.1. Assume that g^ is a complex-valued function and that there are r^ s > 0,(5 G 
(0, tt/ 2) such that, 

' n — S 


9h 


5 

h 

-S 


9h[j^+iy 


ly 


\9h{k)\ < Clk 


9h 

9h 

1 

sK^ 


h 

TT + (5 


+ iy 


h 

Vfc e z. 


ly 


^Ce^sinSsinhiyh)^ /or all y < 0 
< for ally>Q 


Then: 

(1) rs < 1 There are nonzero functions that satisfy the hypotheses. 

(2) rs = l^gh = (and gh{k) = Ch{ir/h^)). 

(3) rs > 1 ^ 5 ,* = 0. 


This is the version we are going to use in order to have a discrete version of Hardy’s theorem. 
As we have pointed out in Section 3, it is not reasonable to think that the condition on the 
coefficients a and (3 that measure the decay of the solution at each time will be the same as in 
the classical result, so we can only talk about weaker versions in this setting. More precisely, we 
have the following result: 

Theorem 4.1. Let fm{t) be a solution to the discrete Schrodinger equation ([ 2 ]), and assume that 
l/fe (0)1 < ) l/fe (1)1 < with a and (3 positive numbers satisfying a + /3 < 2, then 

= if^) i 0. 


Proof. We define the function gh{x, t) such that gh{k, t) = fl({t). As we have seen in Section 2, 
the decay of the initial datum implies that the periodic function gh can be extended to a entire 
function and that 


\9hiz,0)\ < 


_ _ _cosh(y/i) 

V^Io^a/h"^) 


for z = x-\-iy. Therefore, solving the discrete Schrodinger equation in this periodic setting, since 
it is a multiplication by another periodic and entire function, we have that at time t = 1 the 
corresponding function gh{z, 1) is going to be periodic and entire as well. Moreover, 


\9h{z,l)\ 


\eh^ 


• cos{zh) 


9hiz,0) \ < 


h 

V^Io{a/h‘^) 


2 

eh^ 


sin{xh) sinh{yh)-\--^ cosh(y/i) 


The key point here is how to choose S in the corollary in order to hold all the hypotheses. 
Since a + /3 < 2, there is e such that a + /3<2 — e, soa<2 — e. For this e we take close but 
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less than 7r/2 such that . = 1 + —. and = 2 -= 2 — a — e>0. Now, if w < 0, 

' Sin dg Q! ^ sin Og ’ ^ — ' 


gh{-^+iyA 


< 


-eh 


h a .yh 


sin (5g sinh(y/i) + ’j^ cosh.{yh) 


y/^Io{ajh'^) 


+ sin<5e sinh( 2 /h) 


< 


v^/nfa/h^) 


+ ^ sinh(y/i) 


SO we have exactly the kind of bounds we use as hypothesis in Corollary 4.1. In order to get 
the result, we may let the constant C in the corollary depending on h and a. However, using 
the asymptotic expression of Io{t) we get rid of the dependence in h, since we have, at least for 
h small enough 

h a 3y/a 

— r= -^ -. 


It is easy to check that for the other three lines in the plane we need to consider we can do 
the same, so gh{z^ 1) satisfies the hypotheses of the corollary with r = r^. On the other hand, 
the decay of the Fourier coefficients at time 1 implies that s = ^ and again the constan C in the 
statement of the corollary is independent of h. Finally, we have that. 


rs 



> 1 


o; + /3 < 2 — e. 


so, by the corollary gh{z, 1) = 0 gh{z, t) = 0. 


□ 


Remark 4.1. The condition on a and /3 is, as it should be expected, weaker than the original 
condition a/3 < 4. Actually a + P<2^ a/3 <1. Going back to the example we considered 
in Section 3 to show that there are nonzero functions satisfying the hypothesis of Theorem 3.1, 
recall that then a = 1/2 and [3 = 5/2, so in that case a + /3 = 3, hence, that example does not 
make any contradiction with this result. 


Remark 4.2. We can easily see that the condition a + /3 < 2 is optimal, in the sense that we 
can give examples of solutions satisfying the hypothesis when a + P = 2. Indeed, consider 
phfa\ ^ 


Jk\^)— 7o(l//i2) 


fHt) = 






^ /„(-^//^2)4_„(2^^//^2) 


/.(l) = e 




/0(l//l2)- 


Hence, we have a non-zero function that satisfies the decay conditions with a = P = 1. On 
the other hand, looking at Theorem 3.1 we have that if the function uq of the hypothesis (0) 
exists, it has to be zero, and it is as easy to check as in the example in Section 3 that this is the 
case. 


Roughly speaking, what this theorem is saying is that in order to have a nonzero solution, if 
one of the coefficients is small, then the best possible coefficient one can have at time 1 is close 
to 2, thinking about being "best" as being the smallest possible coefficient, since the smaller the 
coefficient a or P, the better the decay. On the other hand what the classical Hardy’s principle 
states is that the smaller one coefficient is, the bigger the other coefficient has to be in order to 
have a nonzero function. 
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Considering now the heat equation, we can use the same argument as for the Schrodinger 
equation in order to have a very similar result, again asking about decay conditions at two 
different times. Now solving the heat equation is equivalent to multiply in the periodic setting 
by so we use Theorem 2.1 with 6 = 0 and 0 = 0 in order to have: 


Corollary 4.2. Assume that is a complex-valued function and that there are r^ s > 0,(5 G 
(0,7r/2) such that, for all y G M, 


9h 


r/2 + (5 


ly 


9h 


r/2 + (5 


ly 




\ghik)\ < ( — ) , Vfc G Z. 


Then: 

(1) rs < 1 There are nonzero functions that satisfy the hypotheses. 

(2) rs = I ^ gh = (and gh{k) = Chir/K^)). 

(3) rs > 1 ^ g?, = 0. 


And by means of this Corollary, we have the result: 


Theorem 4.2. Let v!f,{t) be a solution to the discrete heat equation ([31), and assume that 
kfc(0)| < I kfe(l)l < ^ ^ positive numbers satisfying a + /5 < 2, 

then = (u^) = 0. 


Proof. The proof of this theorem follows the same argument as Theorem 4.1, just rewriting it 
for solutions to the discrete heat equation instead of the discrete Schrodinger equation. 


□ 


Moreover, we can see again that the condition o; + < 2 is optimal, but now we can construct 

examples that are not discrete versions of the identically zero function and are close to the 
optimal condition a P = 2. For instance, if we take the function uo{x) = e~^ then 
u{x,t) = e^^^^uo{x) is the function given by 


u(x, t) 


g-a;V(4t+2e) 

^Jltf e 1 


so in this case, we can take 5 on Hardy’s uncertainty principle for the heat equation bigger than, 
but as close as we want to 2 + e. Then, if we take a discrete version by using the modified Bessel 
functions, we have a result with a = e and /3 = 2 + e, so o: + /3 is as close as we want to the 
optimal value 2. 


5. Appendix: Examples 

In this Appendix we are going to see some examples of nonzero functions that satisfy the 
hypotheses of Theorem 2.1 when we are in the first case of the theorem rs < 1. 


First example 
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The first example, and the most basic example, is given by the function fh{z) = Ce'^ cos{zh b) ^ 
where v = oe*^, for r < a < 1/s. 

Second example 

A calculation shows that the Fourier coefficient of Ce'^ is Ik{u/h'^)e~'^^^. Here we 

are going to take u such that |m| = 1 and a > 1, in order to define the function whose Fourier 
coefficients are 

Mk) = (ae--ri. (i) . 

This function satisfies the hypotheses for s = 1 and u = Indeed, one can compute the 
function fh to get 

and now it is easy to check that all the hypotheses are satisfied. 

More examples 

In some particular cases, like the particular cases of Corollary 4.1 or Corollary 4.2, we can get 
more examples by studying the evolution of discrete equations. For instance, if we assume that 
we are in the case of the discrete heat equation, that is, we have b = 6 = 0, and, for the sake 
of simplicity, we assume that h = 1, once we have an example of a function that satisfies the 
hypotheses, say /, we can generate more functions of this type by using the original function as 
initial datum of the discrete heat equation. Indeed, for the evolution of the Fourier coefficients 
we have 

f {k,t) = ^ /(m,0)/m_fe(2t) ^ \ f{k,t)\ < 4 + 2t 


On the other hand 

\f{-K/2 + 5 + iy,t)\ = e-2*-(2i+Osin5cosh(y)^ 

SO the function f{z, t) also satisfy the hypotheses of the corollary, with rt = r + 2t and st = 

Notice that since rs < I, we have immediately that rtSt < 1 for all t. 

Furthermore, if we have a function not necessarily on the class of functions that satisfy the 
hypotheses, but it satisfies the condition on the Fourier coefficients, let us say a function g such 
that for some a > 0, 

|ff(fc)| < Ch{a), 

and we solve again the discrete heat equation then we have that 

\g{k,t)\ < CIk{t + a), 

\git,x + iy)\ < 


If t < a we do not have the existence of 6 such that the assumptions hold. However, if 
t > 2a, there exists 6 such that t cos(7r/2 + 6) = a(cos(7r/2 + 5) — 1) (Notice that cos(7r/2 + 6) = 
cos(37r/2 + d)). 

Hence for those t we have that 

\fit,Tr/2 + 6 + iy)\ < coshy^ 

So f{t) is in the class of functions which satisfy the assumptions of the corollary. Notice that 
the product of the two constants is < I. 
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On the other hand, we can fix S and then, when t > a ^ the hypothesis holds. 
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